Excitations of Bose-Einstein condensates in a one-dimensional periodic potential 
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We report on the experimental investigation of the response of a three-dimensional Bose-Einstein 
condensate (BEC) in the presence of a one-dimensional (ID) optical lattice. By means of Bragg 
spectroscopy we probe the band structure of the excitation spectrum in the presence of the periodic 
potential. We selectively induce elementary excitations of the BEC choosing the transferred mo- 
mentum and we observe different resonances in the energy transfer, corresponding to the transitions 
to different bands. The frequency, the width and the strength of these resonances are investigated 
as a function of the amplitude of the ID optical lattice. 

PACS numbers: 67.85.Hj, 67.85.De 



The knowledge of the linear response of a complex 
system gives crucial information about its many-body 
behavior. For example, the superfluid properties of a 
three-dimensional (3D) Bose-Einstein condensate (BEC) 
are related to the linear part of the phononic dispersion 
relation at low momenta [l|. The presence of optical 
lattices enriches the excitation spectrum of a BEC in a 
remarkable way. For deep three-dimensional lattices, the 
gas enters the strongly correlated Mott insulator phase 
and the spectrum exhibits a gap at low energies Q • The 
response of a BEC in the superfluid phase is also drasti- 
cally modified by thepresence of a one-dimensional (ID) 
optical lattice % H i, i, Q . Indeed, as in any periodic 
system, energy gaps open in the spectrum at the multi- 
ples of the lattice momentum and it is possible to excite 
several states corresponding to different energy bands at 
a given value of the momentum transfer [H, In addi- 
tion, the linear dispersion relation of the superfluid, and 
thus its sound velocity, is changed. In the mean-field 
regime of interactions these peculiar features of the exci- 
tations of a superfluid BEC in the presence of an optical 
lattice are captured by the Bogoliubov theory [l[ . 

Bragg spectroscopy represents an excellent experimen- 
tal tool to investigate the linear response of gaseous BECs 
[HI ] . It has allowed to measure the dispersion relation of 
interacting BECs in the mean-field regime [111, [IH, EH , to 
characterize the presence of phase fluctuations in elon- 
gated BECs 14 1, to study signatures of vortices [lj| and 
to study strongly interacting 3D Bose [l6| 



14j . to study signatures of vortices [If 

more recently 

and Fermi [17[ gases close to Feshbach resonances as well 
as ID Bose gases across the superfluid to Mott insulator 
transition |18j . 

In this work we use Bragg spectroscopy to probe the 
excitation spectrum of a 3D BEC loaded in a ID optical 
lattice. Previous experimental studies have so far in- 
vestigated the excitations of superfluid BECs within the 
lowest energy band of a 3D optical lattice by means of 
lattice modulation [l!| and Bragg spectroscopy [H, H(| • 
This paper presents a detailed experimental study of the 
different bands in the excitation spectrum of an interact- 
ing 3D BEC in the presence of a ID optical lattice. We 
measure the resonance frequencies, the strengths and the 
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FIG. 1: (a) Measured BEC excitation spectrum in the pres- 
ence of a lattice with height s = (22 ± 2). The increase Act of 
the width of the atomic density distribution is monitored as a 
function of the relative detuning v between the two counter- 
propagating Bragg beams. The data are fitted with Gaussian 
functions (black line). The arrows below the resonances in- 
dicate the corresponding bands, represented in (b) with the 
same colors, (b) Band structure of the excitation spectrum of 
a BEC in a ID optical lattice with s — 22: first, second, third 
and fourth bands are represented (gray, green, blue and red 
lines respectively) . The arrows indicate the processes starting 
from a BEC at q = and inducing the creation of excitations 
in the different bands at a quasi-momentum transfer 0.12g_L. 



widths of the transitions to different bands of the ID op- 
tical lattice. The measurements are quantitatively com- 
pared with Bogoliubov mean-field calculations for our ex- 
perimental system 

We produce a 3D cigar-shaped BEC of N~ 3 x 10 5 
87 Rb atoms in a Ioffe-Pritchard magnetic trap whose 



2 



axial and radial frequencies are lo v = 2n x 8.9 Hz and 
lu x = lu z = 2t: x 90 Hz respectively, corresponding to a 
chemical potential /x ~ /ixl kHz, with ft, being the Planck 
constant. The condensate is loaded in an optical lattice 
along the longitudinal direction (ij axis). Two counter- 
propagating laser beams with wavelength Al = 830 nm 
create the lattice potential V(y) = sEr sin 2 (q^y) where 
gL = 27t/Al = 7.57 iim -1 is the wave-number of the 
beams, s measures the height of the lattice in units of 
the recoil energy Er = h 2 /2m\'[ ~ h x 3.3 kHz, m is the 
mass of a 87 Rb atom. The loading of the BEC in the lat- 
tice is performed by slowly increasing the laser intensity 
up to a height s with a 140 ms long exponential ramp 
with time constant r = 30 ms. 

After a holding time of typically 20 ms in the lattice 
at the height s, we excite the gas by shining two off- 
resonant laser beams (Bragg beams) for a time Ats = 3 
ms. The Bragg beams induce a two-photon transition 
transferring momentum and energy to the atomic sam- 
ple. Their wavelength is Xb = 780 nm corresponding 
to a wave-number qs = 2tt/Xb = 8.05 /Ltm , and they 
are typically detuned by 350 GHz with respect to the 
D2 transition of 87 Rb. To change the transferred mo- 
mentum we use two different geometries of the Bragg 
beams. In the first configuration the two beams are 
counter-propagating along the y direction and the trans- 
ferred momentum is q = 2qB = 2.\2q^, that corresponds 
to a quasi- momentum 0.12g£. In the second configura- 
tion the angle between the Bragg beams is smaller and 
the measured value of the transferred momentum (and 
quasi-momentum, in this case) along the y direction is 
q = 0.96<7l. In both the cases, the two beams are de- 
tuned from each other by a frequency difference v using 
two phase- locked acousto-optic modulators. We quantify 
the response to the excitation by measuring the energy 
transferred to the gaseous BEC. The measurement of the 
energy transfer E(v, q) is connected with the dynamical 
structure factor S(u, q) (giving information on the exci- 
tation spectrum) by the relation [l[ 

E(v,q) <x vS(v,q), (1) 

where v and q are the frequency and the momentum of 
the excitation. In particular, this result applies for long 
enough Bragg pulses, namely v/S±b ~> 1, which is the 
case in our experiment since v is of the order of several 
kHz. 

In order to get an estimate of the transferred energy 
E(v,q), we adopt the following procedure. We linearly 
ramp in 15 ms the longitudinal optical lattice from s (the 
lattice height at which we have applied the Bragg pulse) 
to the fixed value st = 5. Then we let the excitation 
being redistributed over the entire system by means of 
the inter-atomic collisions for 5 ms. After this interval 
time we abruptly switch off both the optical lattice and 
the magnetic trap, letting the cloud expand for a time- 
of-flight ixoF = 20 ms and we then take an absorption 
image of the density distribution integrated along the x 
axis. Since the atoms are released from an optical lat- 



tice of relatively small amplitude (s/ = 5), the density 
distribution exhibits an interference pattern [2l[ . We ex- 
tract the rms width ct of the central peak of this density 
distribution by fitting it with a Gaussian function. The 
increase Act of this quantity is used as a measurement of 
the energy transfer [19( . For a given value of the trans- 
ferred momentum q and amplitude s of the lattice, this 
procedure is repeated varying the energy hv of the exci- 
tation in order to obtain the spectrum. In our regime of 
weak inter-atomic interactions, the excitation spectrum 
of the BEC in the presence of a ID optical lattice can be 
described by the mean-field Bogoliubov approach [U, [f| , 
by which wc calculate the resonance frequencies Vj and 
the transition strengths Z- } to create an excitation in the 
Bogoliubov band j. 
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FIG. 2: (a) Band spectroscopy of a BEC in the presence of 
a ID optical lattice: the energy of the resonances is reported 
as a function of the height s of the lattice. The experimen- 
tal points (green squares, blue circles and red diamonds) are 
compared with the numerical calculation of the Bogoliubov 
spectrum in the presence of a ID lattice (solid lines) and the 
single-particle Bloch spectrum (dashed lines). The lines cor- 
respond to the energy of an excitation in the second (green 
line), third (blue line) and fourth (red line) Bogoliubov bands. 
In the inset of (a): zoom of the graph (a) for low values of 
s. (b) Relative strength of the excitations in the j = 2, j = 3 
and j = 4 bands. Symbols and colors are the same as in (a). 

We first discuss the results obtained with the configura- 
tion of counter-propagating beams, i.e. for a transferred 
momentum q = 2.\2q^. The induced two-photon tran- 
sition is characterized by a measured Rabi frequency for 
the BEC in the absence of the optical lattice ~ 2-7T x 1 
kHz for the typical power and detuning of the beams used 
in the experiment. A typical Bragg spectrum is presented 
in Fig.QJa) corresponding to a lattice height s = (22 ±2). 



The spectrum exhibits multiple resonances corresponding 
to the creation of excitations in the different Bogoliubov 
bands as shown in Fig. QJb). From Gaussian fit of each 
resonance we extract the central frequency, the width 
and the relative strength of the transition towards the 
corresponding band. In Fig. EJa) , we plot the energy val- 
ues corresponding to the measured central frequencies as 
a function of s. The vertical error bars come from the 
result of the fitting procedure while the horizontal error 
bars correspond to possible systematic errors in the lat- 
tice calibration (estimated within 10%). For large enough 
amplitude s of the periodic potential we observe up to 
three different bands. We find a good agreement between 
the experimental data and the numerical results of the 
Bogoliubov calculation (solid lines in Fig. [21(a)). In par- 
ticular, for low amplitudes of the ID lattice (s < 6) the 
agreement of the resonance energies with the Bogoliubov 
bands (full lines) is better than with the single-particle 
(dashed lines) Bloch bands (see inset in Fig. [21(a)). For 
larger amplitude of the ID lattice, we can not explic- 
itly distinguish between the Bogoliubov and Bloch re- 
sults. This comes from the experimental uncertainty on 
the calibration of the lattice amplitude. 

In the entire range of s values used in this work, we ob- 
serve a resonance corresponding to an excitation created 
in the third band j = 3 (blue circles on Fig. For larger 
lattice amplitudes two other resonances appear, respec- 
tively for s > 4 and s > 20, corresponding to an excita- 
tion in the j = 2 band (green squared in Fig.[2ja)) and in 
the j = 4 band (red diamonds in Fig.[^a)). This demon- 
strates the possibility to excite, in a periodic system, sev- 
eral states for a given momentum transfer 8]. For weak 
optical lattices creating an excitation in the j = 3 band 
is the most efficient process since the excitation energy 
of this band is continuously connected as s — > to that 
of the BEC in the absence of the ID optical lattice at 
the transferred momentum q — I.Ylq^. On the contrary, 
the possibility to excite states in the second and fourth 
bands of the optical lattice requires a large enough am- 
plitude s. These observations can be quantified in terms 
of the strength Zj of the different excitations, which can 
be extracted from the energy spectrum. The strengths Zj 
are proportional to the integral J dv Sj (q, v) with Sj (q, v) 
being the structure factor corresponding to the creation 
of an excitation in the Bogoliubov band j [j| . From Eq. [1] 
and assuming that Vj is much larger than the width of 
the resonances of Sj(q, is), we obtain 

z j(<7) J dv Sj{q,is) oc ^- J dv Ej(q,u) = g- } . (2) 

In the experiment, we extract the quantity gj from a 
Gaussian fit of the different resonances. Normalizing the 
sum of these quantities to one for the first three observed 
resonances (<?2 + .93 + .94 = 1) allows direct comparison 
with the relative strengths ZjjiZ^ + Z$ + Z4) for j = 
2,3,4. The comparison between the experimental data 
and the calculation reveals a reasonable agreement (see 
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Fig. [Kb)). 
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FIG. 3: Rms width of the resonances j = 2 (green squared) 
and j = 3 (blue circles) as a function of s. The gray region 
corresponds to the experimental rms width (with its uncer- 
tainty) for the BEC in the absence of the lattice (s = 0). 
The green and blue lines are respectively the bandwidth of 
the second and the third band, calculated in the mean-field 
Bogoliubov approach. 

From the Gaussian fit of the experimental spectra (see 
Fig. [He)), we also extract the rms width of the reso- 
nances j = 2 and j = 3 with the results plotted in Fig. [3] 
Different sources contribute to broaden the observed res- 
onances. The inhomogeneous density of the trapped BEC 
is a first source [11| . From the measured spectrum of the 
BEC in the absence of the optical lattice (s = 0) we 
extract this contribution being (0.36 ± 0.11) kHz, con- 
sistent with the expected value ~ 0.26 kHz JTTJ - The 
other contributions to the width are related to the Bragg 
spectroscopic scheme. For our experimental parameters 
the largest contribution comes from the power broaden- 
ing (Aisp ~ 1 kHz), whereas the atom-light interaction 
time broadening (Ais t ~ 167 Hz) is much smaller. The 
total resonance width can be obtained by quadratically 
adding up all these rms contributions. In the presence 
of the optical lattice we observe that the widths of the 
resonances corresponding to the excitations in the bands 
j = 2 and j — 3 lie within the experimental range of 
the resolution as expected for a coherent system, except 
in the case of j = 3 for large amplitudes of the lattice 
(s > 20) where the width is much larger. We attribute 
these larger widths at high amplitude of the ID lattice 
to the long tunneling times (0.11 s for s = 20) implying 
that the system is not fully coherent along the y direc- 
tion on the time scale of the experiment. Indeed, the loss 
of coherence spreads the population of quasi-momcnta 
across a larger fraction of the Brillouin zone. This re- 
sults in a wider range of resonant energies in the system 
and, for large amplitude s, one expects the width of the 
resonances in the energy spectrum to be equal to the 
bandwidths. In Fig. Owe have plotted the bandwidths 
of the j = 2 and j = 3 bands (green and blue lines). 
When the system becomes incoherent the width of the 
resonance j = 3 is equal to the bandwidth. This effect is 
not observable for the band j = 2 where the bandwidth 
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(green line in Fig. [3]) is smaller than the experimental 
resolution. 
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FIG. 4: (a) Excitation spectrum of a BEC in the presence of a 
ID lattice with height s = 11 at a transferred momentum q = 
0.96(71, along the y direction. Fhe arrows below the resonances 
indicate the corresponding bands, represented in (b) with the 
same colors, (b) Energy of the resonances corresponding to 
an excitation in the band j = 1 (gray circle), j = 2 (green 
squared) and j = 3 (blue circle) as a function of the transferred 
quasi-momentum q for a fixed value of the lattice height s = 
(11 ± 1). Fhe experimental points are compared with the 
numerical calculation of the energy bands in the Bogoliubov 
approach for s = 11 (gray, green and blue dotted lines); the 
solid lines correspond to the bands for s = 10 and s = 12 to 
take into account the 10% uncertainty of s. 

We also perform the experiment with a different con- 
figuration of the Bragg beams corresponding to a trans- 
ferred momentum along y q — 0.96<7l. In Fig. IDJa) an 
excitation spectrum in the presence of an optical lattice 



of height s = 11 is depicted. Note that a first resonance 
at low frequency is visible corresponding to an excitation 
with non-zero momentum within the lowest energy band 
(j = 1). Such a resonance is not observed using counter- 
propagating Bragg beams because the strength of this 
transition is negligible for q = 2.12qL. Due to the varia- 
tion of the transferred momentum q with respect to the 
previous case, the frequency of the resonances are shifted 
according to the dispersion relation of the different en- 
ergy bands of the system. In Fig. 0Jb) we report the 
frequency of the resonances j = 1, j = 2 and j = 3 for the 
two values of quasi-momentum used in the experiment 
(0.12c7l and 0.96<7l). We use the region comprised be- 
tween the calculated bands for s — 10 and s = 12 (solid 
lines) to take into account the 10% error in the lattice 
calibration. The experimental points are in good agree- 
ment with the numerical calculation of the Bogoliubov 
bands for s = 11 (dotted lines in the Figure) within the 
experimental uncertainty. 

In conclusion, Bragg spectroscopy has been used to 
probe the response of a Bose-Einstein condensate in the 
presence of a ID optical lattice. Changing the angle of 
the Bragg beams allowed us to investigate excitations for 
a transferred quasi-momentum close to the center and 
to the edge of the reduced Brillouin zone. We have ob- 
served different resonances in the response function of the 
system corresponding to the different bands of the peri- 
odic potential. Being the system in a weakly interact- 
ing regime, the experimental results are in quantitative 
agreement with the Bogoliubov bands approach. This 
work opens the way to investigate the modification of 
the excitation spectrum in the presence of an additional 
lattice with different wavelength (bichromatic potential) 
[22| and eventually to study the localization of the exci- 
tations in a true disordered potential [231 ]. 
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